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ABSTRACT 


Pwo tramsport equation has been developed that 
the bupome dwseribution, aS a function of position, 
VelOcliyyeana size, must satisfy. An analytical model for 
bubble transport in the upper ocean is chosen and solutions 
Pa wdewelopeamtorerniIs model. For a surface source only, 
inese Solutions Compare favorably to experimental data in 
ime leap -suUreace region. “The depth of this region of 
agnecmeme depends upon the circulation field chosen, but 


is of the order of 3 meters. 
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DESCRIPTION 

depth (z=0) 

bubble welocity 
Circulation pieidevelocaty 
bubble radius 

lil yale 


bubblesdmstadbutaon 


(functiom of v,r,1) 
distributive source 
bubble volume 

density of water 

density G2 vas in Dubbie 
viscosity of water 
relative velocity (v-V) 
Dubble acceleration 


fractironat pare of bubble 
volume 


9/2 n/op,, 
atmospheric pressure 
hydros atime pres sum ec 
gravity 


bubble distruputaen 
(Guamction of ZyX5V_sVi 1) 


Po /pwg | 
bubble density (as a func- 


fionoe 5251), | ¥d>v 
position 





I. INTRODUCTION 


Studies concerning the prepagationsof sound in the 
ocean and under water sound scattering from the sea surface 
have shown a discrepancy between existing theories and ex- 
pPerimnemeal results. Because airvhas asWarkedlysditferent 
density and compressibility than sea water and because of 
Micwmcoemant Characteristics of bubbles, the suspended air 
bubbles in sea water have a profownd effect upon underwater 
sound and can explain part of the discrepancies that have 
been encountered (Ref. 1 and 2). There are theories that 
predict the effects of air bubbles on absorption, scatter- 
im@eana atcenuatrzon of sound as well as their effects on 
fommemvelocity (Ref. 3), but there 1s no theory at present 
mimem preaters the distribution of air bubbles in the 
S@eam, ovch a theory would enhance existing theories and 
mom: 1 better predictions on the behavior of sound prop- 
agation in the upper ocean. 

A few experimental studies have been conducted on the 
dastimpution of bubbles in the ocean. Blanchard and Wood- 
POCl@cmceploncers Im this £1e1d, using a bottle to scoop 
up water from breaking waves and counting the number of 
bubbles of eeenent sizes (Ref. 4). Later studies con- 
ducted at the U.S. Naval Postgraduate School, Monterey, used 
leeucelentcchmmques to determine the bubble population, 


taking advantage of the fact resonant acoustical scattering 





and absorption cross section of bubbles are of the order 
of a thousand times the geometrical cross section. These 


experiments indicated that bubbles of radius 60 microns or 


a 
‘3 


@reacer have a distribution depth dependence of Z 7%, 


whereas smaller bubbles follow closer the exponental law 
e 2/L where L is between 5 and 9 meters (Ref. 5). 

In this paper, development of a bubble transport 
equation has been outlined that the bubble CPSieet DUMeOn ads 
Poewmeeron OF pOSIlE10n, velocity and size, must satisfy. 


An analytical model for bubble transport in the upper ocean 


is chosen and solutions are developed for this model. 





II. TRANSPORT EQUATION 


Peter sere approach 1s used here to develop the 
bubble transport equation, since the results are the same 
as With a more rigorous approach and the mathematics is 
Messeinvolyed (Ref. 12). Consider the dastribution 
¥(r,v,1,t)d°*rd°vdl which defines the number of bubbles 
with radii in dl about 1 and with velocity in d>v about v 
which are in a volume d°r about r at time t. If eollisvens 
Mchomenc wuomlmes aremeneclected, then the following equa- 


Pr Olcmecme pec wiatten for the seven-dimenstomal space: 


gt (r,v,1,t)d°rd°vdl = |net fel Ux Slee te ES d>rd5vdl (1) 
where S(#,v,1,t)d°rd°vd1 is a distributed source of bubbles 
introduced in d°r about # at time t with radii in dl about 
1 and velocity in d>v about v. 

Meterring to figtme 1, whach represents a volume ele- 
Mmemeeit real space, the met flux in, per unit volume, for 
mice xeadirection, 1S: 


fluxin-fluxout 


_Vx (x) ¥ (x) dydz-vx (x+dx) ¥(xtdx)dydz a(v,¥) 
Mo Ses ic 


Dinar, the met flux in per Unit volume for the y and 
Z-GieeeGtELOns 15 - and - aC 2) Vespeecimely . 
femeamnine the above results, the total net flux in per unit 
WOMUncCumOn real Space 1S = Wee), where the subscript r 
one the dell-eperator means the operatiom 1S taking place in 


Hrewmeat space. A Similiar approach is used to obtain the 


total net flux in for velocity and size space: - V: (a¥) 











Vy (xt+dx) ¥ (x+dx) 
dydz 
tH» (Flux out) 


Vy (x) ¥ (x) dxdz 


(Flux in) 
AZ 


Fisoune 1 


Elemental Volume Showing "Bubble" 
Streaming in Real Space (x,y,z). 





for velocity space and - d(RY) for size space, where R=dl/dt. 
o1 


@ne total net fluxein for the seven-dimensional space under 
Bonsideratiomeis thessum of the net flux in for real, 
MelGeity, and sizé spaces. Substituting this into equation 


(1) gives the general bubble transport equation: 
ev!) = V_- (aY 9 (RY 
ge. - ¥,. (VP) - V,-(8Y) - 2RY) , i 

miceremaimaer of this paper will consider the special 


eice Or Steady-State with no distributed sources, for which 


equation (2) can be written 





- = Ray _ - oR 
ae V Y 7 e V Y as or, aoa e oS 
- a.¥ a lv a + S| Me oD) 
pace rand v are Miaecpendent Vartables. Although 


distributed sources will be neglected in this thesis, 
ierace SOuUTrceS can be introduced as boundary conditions. 
it should be pointed out that the non-zero right hand 
side of Cipeenon (5) 1S peculiar to the bubble transport 
probiem and arises becouse a is a function of v and R 


Mra tunction of 1. 





III. ANALYTICAL MODEL FOR TRANSPORT EQUATION. 


PeOmemoolutions tO equation (3) can be obtained, a 
Pembiartc model for the bubble transport equation must be 
dewemoped, with particular attention to formulating ex- 
pressions for the bubble acceleration (a), the time rate of 
Grange sot the bubble radius {(R) and the circulation field 


(Vv). 


A. MODEL ASSUMPTIONS 

Mie sree lMused in this development assumes, Reynold's 
fompers Less “than one, which is the Stokes' Law Regime. 
With regard to this assumption two comments should be made. 
fparsct, tor small Reynold's numbers, bubbles will behave 
Misemsclid spheres with little or no distortion and secondly, 
wes places am Upper limit of about 150 microns on the size 
Ciumembpbles that cam be comsidened. For Reynold's numbers 
iuemecceater than one, (larger bubbles) the resistance of 
epewturbulent wake will have to be included in the equation 
Cer nron, alone with Stokes! drag force (Ref. 6). 

Surface tension has been neglected in the model, so the 
pressure inside the bubble is equal to the hydrostatic 
mibtdepuessure. since the pressure from surface tension 
2y¥/1, where y is approximately 74x10°5 NT/m for sea water, 
ms aastimpenon will introduce an error of less than § 


Dem@enmtatoOr radii greater than 30 microns. Thus, where the 


10 





Reynold's number places an upper bound on the radius, neg- 
lecting surface tension places a lower bound on the bubble 
Paaiussunder consideration. 

The volume of a gas bubble in a liquid varies primarily 
HiGomonethewerLeécteote@tiwo influences. Firstj'gas contin- 
uously leaks out of the bubble to dissolve in the surround- 
ing liquid. If the bubble was stationary, the ensuing loss 
of mass would necessarily lead to its ultimate disappearance. 
Onieene other hand, as a bubble assends or desends in the 
ieaguaaed , ssheeanbacnt hydrositatr»e pressure decreases or in- 
Greases! and so doesethe pressure within the gas bubble, thus 
Sereceting thems ize% 

iitcmeonhoe ts of SeReminteractron= between gasedirtiusion 
enuie changes inehydnositta tic spresisure can best be seen by 


looking at the ideal gas law: 


an13 = rITN/P 


mirere N 15 the total number of gas molecules, r is the gas 
Soitemiioe | 15 the temperature, and P is the pressure. Now, 
balamme the time derivative of the id@al gas law equation 


and assuming isothermal changes, 


2 l 
4n1°R = . PdN/dt - Ndp/dt 


amc Lrearwangaine : 


1 dN/dt dP/dt 
3 N Ie 


Jal 





The above development assumes a single gas inside the 
bubble vice some mixture of gases. In eqmation (4), the 
first term is a result of gas diffusion whereas the second 
MP ieomimeswin une eftectss of compress vem 

mes tO the relative motion between tle water and bubble, 
assume that the composition of the water in respect to dis- 
solved gas is everywhere uniform and saturated at atmospheric 


Pimeosure, P Exeepu in a thin shell surreemding the bubble. 


0? 
Assume that the gas in the bubble is uniform right up to 

themeas water interface, and that the liquid in contact 

focheehe bubbleseis saturated at the pressure, P, inside the 
bubble. Further assume that diffusion gradient is uniform 
Pimeuchout the thinsshell. With these»assumptions, Fick's 
First Law can be used, which relates the amount of material 


which diffuses through an unit area to the product of the 


Gomecmeration gradient and diffusion constant: 


dNdt = - 64n14(P-P_). 


In this expression 6 is a proportionality constant which 
@emencs son d, the thiekness of the shell; on A, the dit- 
fusion constant of the gas in the water; and on a, the solu- 
Dmiitiy eo une oas in the water, in accordance with the 


wewace ron : 


6 = Aa/d. 


eZ 





Weamoervek™s First Law with the hydrostatic pressure rela- 


tionship , Pep. - 0,82; equation (4) becomes: 
j Vv 
ea. ee (4) 
(aj-z) 3 (a, -z) 


where o1,= P,/pyg and v, is the bubble velocity in the 
Boone ction. in neference 7, the expression érT was deter- 
Mimedeexperimentally to be approximately a constant, 


10~° m/sec. 


B. EQUATION OF MOTION 

Eiewme air bubbles are entrained im a Circulation field, 
there are three forces acting on them along with the bouyant 
pemec.  hiscure 2 illustrates am air bubble in a circulation 
field and shows graphicly the various forces which need 
to be considered in order to determine the bubble 
aeccleration.. 

Pietcicmoupmbes were compbetely entrained by the fluid, 
ies, bUbDIC yvelloe1tye= Cimcullation field velocity, they 
would experience the same force as would the fluid enclosed 
in the same volume, p,vdV/dt. However, complete entrain- 
ment Can only occur if o = p,, which is far from being true 
for air bubbles. Since the bubbles are only partially en- 
traaned, the fluid flows past the bubbpites Creating a rela- 
tive velocity (u) and thus introduces a drag force. As 
mentioned earlier, small Reynold's numbers have been assumed 
and thus the drag acting on the bubble can be feeomned by 


Stokes' Law, 


igs 





Bouyant Force 
‘ Drag Force 


A 
,Adhering Mass 
Force 


Entma@inment Force 


Payoumes 2 


Free Body Diagram of a Bubble, 


Shewine thesrorces Acting on 
ties Bubble an a circwlat ion’ -F ive ld 
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Drag = Gaim. 


tiemehirdy torce acting on the bubble is the innertial 
adhering mass force. When a spherical shape moves through 
the water, it drags along an amount of water equal to a 
ireeimonal amount of the volume of the Sphere. In the 
Govations that follows, sigma (co) represents this fractional 
Pavey eameeein tererence (8), sigma has been derived as 0.5 
Lon tansphere. 

Referring to figure 2 and using Newton's Second Law of 
fimo eene Cauation of motion for™an air bubble in a cir- 


Gumaitiom’ field can be written as: 


pv gv. = ——— Dilogse = Og esc (po, -p)¥gk. 


Using p<<p,, and dv/dt = (v'V_)V yields 


a =F Wee Bk + AoW Vr) i Cy 
dl 


where a = 9/2 n/ 9p, and u=v-V. 


eo CRRCULATION FIELD 

PMuamwiomise known about the circulation £ield in the 
upper ocean. Although some work is being conducted at pre- 
Bememac che U.s. Naval Postgraduate School, Monterey. In 
mememenecm( Oo), Levich develops a two-damensiomal equation 
Which describes the particle motion by wave action in an 
mae werleciid=e sand 1S the basic form used here to describe 


erreulation field: 


15 





ae Voie <?sin kxi - Ve *X2cos kxk. (6) 
irs E€QUation gives a periodic form for the horizontal 
direction as well as an exponental decay in depth, which 
might be expected Home tne piysical Sitiation in) the upper 
ecm. I(t 1s also well to note that the condition for an 
imeempressible fluid has been satisfied with this equation, 
meow ve V—0. 

Froure 3 shows the cell structure implied by this cir- 
culation field. For this model, any two adjacent cells are 
Picitimedintewevery Other pair of adja@ecient cells that are 
fomamaeed from Che first pair by an integer number of 
wavelengths. This implies that the boundary conditions must 


be periodic unless unphysical surface sources are to be 


assumed. 
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reguie .3 


C21 Structure of the’ Circulation 
Field Model 


lee 





Tv) SOLUTION TO TRANSPORT EQUATION 


Pieeitomece tion, a look 1s taken at the bubble transport 
equation with the model that was presented in the previous 
section, The main interest will focus at determining the 
bubble distribution as a function of depth and bubble size. 

Diese ince tle bubble transport equation the 
expressions £or a and R from the preceeding section results 


in the following equation: 


viVp¥ + a-Vy¥ + Rav/al = [3e/1 : ¥2/3(a, ~2)| : 


Ey adetining 


[ee 


$(2,X,V,.Vyo1) = |] ¥CE.V,1)ayavy 


OO 


and using the two-dimensional circulation field, the 


above equation can be integrated over the y and Vy Gasre citaion Ss 


and the problem reduces to solving the partial differential 


equation: 


v.99 Vater’ a, o> ad? Rodd _ {2 x 
x =" oe eat + |e fey] ee 


with the boundary condition $(0,x,v SU on 


1,) corresponding to a surface source. 


x? 
Vxo? 


P@ueerone(7 leis a first order, linear partial differen- 
tial equation in 5 variables and can be solved by the method 
of characteristics (Ref. 9). Using this method, the partial 
differential EqUatelonemssecaquivalent to thie tollowine S char- 
Meeeimst ic Ordinary differential equations in terms of the 


parameter z, which was chosen for mathematical convience: 
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Oriedecerastic Diff. Eqns, Characteristics 


Z<0 Z=0 
ax ‘x - 
dz Vo oe 
dv ax . 
dz 7 Vv_ Be XS 

Z 

dv7= az S 
ae Vy om 26 
dl _R 1=1, 
az. Vy 
dq 20 " 
Hee = | aa - = ¢ = 5 ee eee Ae 
dz ra haeah = $0 ( Or EO. ZO 6) 


A close look at the characteristic ordinary differential 
EG@uatwons Teveals that they ane simply the bubble dymamic 
equations. 

porvame these equations for ¢ would be straight forward, 
Ewecot the equation for Vie is mon-linear and analytic solu- 
mienscenave been sound only for very special cases (Ref. 10). 
nis, approximate me@hods mustebe consideredyto obtain ue 

The method used for this paper to solve for Vv, was to 

Z 


assume Vv = oe us Appendix A outlines this method with 


@ine SOlUtION for U_ as: 


lr" — dz' | [79 | v7 
uy =e /o l’v, of&Jo 7 dz" fg 4 Li_lx 8V7 av, dz'+u 
ve OV. O\V, Ox aeaZe ZO 


Perec clemdt ae acceleration on the bubble is large, the 
bubbles will be entrained by the fluid within the first few 


centimeters of travel and u, thus approaches a terminal 


IES) 





[emoecrty., hererring to the computer output for Us shows 


tiat Che -assumption Ne >> u, can be made and also ap 1s 


Z 
independent of es thus implying that we, MS Didie pe mare t 
of Neve: 
Byway simiaiivar approach, it 1s shown in Appendix A that 
Vy, = Vy. Now, the characteristic equation involving x can 
be written as: 


Gixees ee 
ee Vices 


aamuisime the chosen circulation field, 

x ~ 1/k sin7’+(e7k@sinkx,). 
iis Cquarron amplies that the bubbles follow trajectories 
tiene Liuid and form streamlines as was indicated on figure 
3. Thus choosing a X, one can deSemuene tie pawn nine het ne 
bubbles follow. 


TMiemsolutiaon of the charaeteristic equation for 1 is 


straight forward with 


(Zz 


1=(1-z/a,)*/ h P Srey 2) Cala eeu 


a 
litessecond term in this €quation 1s a result of gas diffu- 
sion and limits the bubble lifetime and hence the depth 
Made SlUrtace @encrated bubbles can penetrate. 
Erommene Nast Characteristic equation, 6 takes the 


general form: 


ie 3 
GU Veg 5720) = 99 (Xos¥yo>Vz9l6) (a-2/a4) EXP 


("20 az 
Oj2y, 


20 





whememd asta function of (X%,Z,Vv vee) 1s obtained from the 


> as 


relationships between the characteristics and their variables. 


Assume a surface source 9, that is separable, 


$5 (Vo Vag to) = XA xo) ¥(VpgQ) Eg); 


andamer convienee,, let 


ac cr ee. Yo) = $(v,,7¥ ). 


xo) Z 


Since v, andmaece, are Weanly independent of v 


3 and v this 


XO 2.0.7 


eieace Will not Strongly affect the solution. Therefore, 


1 
Oe yz,V V1) = X LG 6tv, V6, -v,) la-z/ay) ey, 


[2a ale! 
1 v 
O V 


Z 


New, sdetine the "bubbie density’ as aefunetwon of pesition 


and size, 


@wrzZ,1) = ([ CH 7a. V Cae lidar elev: 
xe oes Saenz 


and integrate the expression for $ over vy and v,: 


= t 7 7 is 
Sieeyzs1) = || Rea ies veer) | (1-2/0) EXP 


220 1 Galea aetna 
GZ 
[obo] x 
To take advantage of the delta-functions, the variables of 
integration are changed, dv dv, = clog. aby 


yy) EROS cl) 1S 


tne? Jocobian: 


an (WoeVer) 


d(v 


J = 
x0?" zo) 


(53 





The eauation of interest now has the form: 


®(x,2z,1) = X(x_JL(1_) (1-2/a,)1/3 Exp | “2 dz'J (8) 
O O Oj} Vz 


Because of the number of mathematical steps involved in 
deriving the Jocobian, the development is found in Appendix 


Mamie tinal approximate result is: 








Wher Cx, 


2 Zia. a ért * 
= + V a 
B oe c- 2) OF 


Thus, the equation for the distribution of bubbles with one 


initial surface velocity can be written as (Ref. 12): 


Danze). = X (Xp) L(19) (1-2/0) 1/5522 ef gaz! (3) 
Peeking at equation (9), the term eP edz Wilieoe che pre - 
dominate feature in solving for ®. Two things should be 
noted about €, First, § 1S seen to be proportional to 1 
when the approximation WN is used. This resvuit illustrates 
tie ampomwtance of selecting a good circulation field model. 
secondly, €asRT and shows that gas diffusions is very 


mportant im the solution for the bubble distribution. 


Zl 





Ve REGUETS 


A eOMputer program Was ewritten for calculating equation 
epee ce G@etcrmime 1f the results were phyeacally realistic. 
In this program, the assumptions were made that Ric ql and 
L(1,)#k1,°7/2 (Kero!) om tice batter assumption, K=T 
was used in the computions since the interest was in the 
omer or the distribution curve rather tivan the actual den- 
meemces. Galleulations for three circulativon field veloci- 


ties were made and the following graphs have been plotted: 


oe Ee as 
Smee) vs. 2 bOr VaEToOuUsS Gadi 1 
Meee vse foriaraousedeptis. 

Z Pcie ax penetracion depen LOr a given Size 


Max 


Sireeewccncrated bubble. All calculations were made for 
X,=0, which implies a straight line trajectory downward. 
ft eseme other X, had been chosen, line integration would 
Mauemce oe perrormea on the characteristic curves in S=space 
(X,25V,.Vz,1). 

mie nimcieed! teciniagque used in compucing the integrals 
was 10-point Gaussian Quadrature, which approximates the 
integrand by a Legendre polynomial series. Although this 
Eeemmndue was tairly accurate for ™most of the integrals, it 
SHmranioe be used in Computing the relative’ velocity “(U7 ) . 
FOE Pomecdmetilations it 1S smoted that the antegrand 15 


ceecmmrally zero except for the last few percent of the 
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Dee ),2) where it increases dramatically. U, was 
therefore evaluated by using double precision 32-point 
Gaussian Quadrature on the latter portion of (0,z). 


MecOmmnetme the relative velocity, the approximation 








Z - 
| cle “5 Bo SE x=0, 
0, ¢v, i Vo 


was used which is accurate provided |Z|< 10 meters. Also, 
Einee the calculations involved taking the difference be- 
miecimevOustien integrals, very little error was introduced. 

It 1S appropriate now to make a few comments about the 
results that have been plotted. Figure 4 is a graph of 
maximun depth versus inital bubble size. This graph shows 
Prreeouibbles Haye a finite life time as a result of gas dif- 
fusion and that smaller bubbles It nOGuced a termmcmsiita Ge 
will not Menetrate as deep “as larcer Bubbles. 

Pepe seo ehnough 7 are plots of © (04221) versus depth 
Home various bubble radii. These graphs show the influence 
mmac tne Circulation field velocity has on how deep the 
Eppes will penetrate. it should also be noted that even 
though the surface density falls off as ae tenet 1S ed 
Sentai) depen where the larger bubbles become more populous 
than the smaller bubbles because of gas diffusion. 

Feeires 8 threugh Il are plots of  (05z,1) versus 
radius for various depths. The main feature about these 
graphs is that the initial distribution appears to persist 
towsome depth which 1s a function of the circulation field 


Velocity . 
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Depth (Meters) 


Radius (Micyons) 





Pcmne 4 


Curves for Maximun Penetration vs... 


Bubble Radius for Various Ci rel 2 tion 
Field Velocities 
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@ dl(number/m*) with dl 





it u 5 


Depth (Meters) 
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Bubble Density vs Depth for Various 
Size Bubbles, Waa 0S m/sec. 
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) dl (number/m°) with dl=l1 micron 


ae 
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1 2 3 4 5 6 
Depth (Meters) 
Fvettee® 0 


Bubble Density vs Depth for Various 
Saagie “Bubb bes , Ue 


=-1.0 m/sec. 
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1 micron 


) dl (number/m°) Wa thedal 


101° 


101° 
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1=404u 


SES ALONG, 


Fomaure 7 


DUMP remDetcwiiywevs. Depthetor Vardous 
Size Bubbles, Vg wee m/sec. 
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® dl(number/m*)with dl 
= 
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Bubble Density vs Bubble Radius 
for Various Depths, v,=-0.5 m/sec. 
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) dl (number/m”) with dl 
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i tees 


1014 





Faeures1 0 


Bubble Density vs Bubble Radius 
eves Wie Ounc e We Pit INS sualieeelecS M/ Sec. 


Sin 





VI. SGONCLUSLONS 


There is little experimental data available to compare 
the results to, however, in reference (5) where measurements 
were taken acoustically, general shapes for bubble distri- 
Pieetoms are grven. One should note that the experimental 
Suede reterence (5) is not for isoradius bubbles since 
the resonant frequency is a function of depth, however, a 
good comparision can be made within a few meters of the 
surface. 

Meeowea the bubble distribution as computed does not 
follow a simple mathematical function, the results for the 
first few meters can be approximated by exponental laws, 
eepenaie om the circulation velocity and bubble size. 

[CCimiCmtimEst atet1eure 5, the distribution drops off 
mmmaly with depth pointing to the fact that for small cir- 
culation velocities, a distributive source is needed since 
Mite esdepth penetration is small. For Vv, 1-0 mysec, the 
Slawes Cal be approximated for the first 353 meters by the ex- 
ponental laws, e7 2/5 for small bubbles and e~2/9 for large 
babiokes. Ihe curves for Vege les m/sec follow fairly close 
e~Z/9 for the smaller bubbles when z2-3.5m and e~Z/15 for 
iieser bUbbLeS wien z2-4.5m. The curves for Voz7i.0 and 


1.5m see compare favorably with the experimental data 


O 
Memenomeds especially for the smaller bubbles. 
Pmenethough there 1s fairly good agreement in the 


distribution near the surface, a definite discrepancy does 


BZ 





Seeitemiemcic deeper depths. Probably, the major reason for 
Merewdiserepancy 1S the @aet that no distributive or bottom 
Pemmeceouwere Considered in the calculatioms. The data in 
reference (5) was taken in shallow water and calm seas, so 
mrelusten Of SUch Sources into the model should improve the 
aeneememe between the solution and the experimental data. 
Another possible reason for the discrepancy is that a large 
portion of the Jacobian was assumed small and neglected. 
iiees term Was not looked at in detail and could be of some 
Maportance at the greater depths. 

Saec see results do detec in part with experimental 
Giese the approach taken in this paper does appear feasible 
Pamcew@ernining the bubble distribution in the upper ocean. 
Tommwnener, Studres in thas areas 1t as suggested that solu- 
tions to the bubble tGiwspereecuuatton bemextemdced to include 
calculations for x#0, models with distributed sources and 
Wemerce tension and@better models for the"eirculation field 


and gas diffusion. 


ae 





APPENDIX A 


Solutions for the Relative Velocities 


In this appendix, solutions are obtained for the rela- 
tive velocities from the characteristic equations. Consider 
Piste the eharacteristic equation for Vo» which is non-linear. 


Write Vics me differentiate: 


paid 
dz dz dz 


PPS teltuting this result into the characteristic equation 


for v, gives 








dz ~*~ 42 vz ~ ovzg + o {vy a, + 5z 


A formal solution to this equation can be written: 


duz, gq Uz g 1 & aVz aVz 


0 7: 
Ove Vz OX OZ 


= 


t 
dz' + 20 ) 


Or wemewrittine, 


iz Zz O dz" 
i = ~),12%z [°° 12y,2 ee | pes ue zs) 
mn e + O 


O 


N 


ol 


[? zZ 1 0. 
Z 


dz" V : 
or ee | B+ ee OVZ 4+ saan a la ae I. 5 az 





OVz O\v, ox OZ 


For computational purposes, the’ approximation v, 2 Wael 
Gemiiem@e Where appropriate, since the above expression is an 
Miplicit equation. This approximation turns out to be valid 


as indicated in the Results Section. 
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iimonder to solye the characteristic equation for Vy > 


the same approach is used as described above, with the 


Peswuit : 
Zz, vagaries 4 
a _.1 Snes | CI OV) 
Ue ed he aoe | a Jo E OV.ox Jae 
We Ive, 1/3V 
° 1(aVx , 
Oo ‘dZ 


V 


2 tafe 


A close look at this equation reveals that u, is approximately 


X 
zero since the exponental term shows a sharp decaying 
feseurewand the circulation field chosen does not c! ige 


rapidly. Thus the approximation V=Vx Gan) De atice de seelie 


Gobet ataons . 
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APPENDIX B 


Dewmelopment of the Jacobian 


irs sectton as devoted to thesmathematical manipula- 
moms involyed in developing an expression for the Jacobian 
momered for €quatvon (8). The folllowang are the functional 
maraevonshiips among the variables and characteristics that 


need to be considered: 


Vy = Vx (2X9 sVxQ0Vz 9219) 
aa AG ot co eee 
1 = 1 (Ae een ll) 
— = xX CARS, VR. » Yo) 
Vx = Vy XQ 0 Vx99Vz0 219» Vx0) 
_ EA OI. <a 


Kor mathematical convience, define: 


Ji = Wey OG 92 = 8Vy/9Vz9 ; Jz = oe Ona 
Jy = dVz/9Vz, 3; J5 = ax/ dies, : Je = 026) Von 
Thus, 
= re Ca) a JiJq a JoJ3- 


Now, using the characteristic equations for Vi, 2V, and x 


and noting that, for example: 


alae) = a ‘ 


Qu 
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fre folrOwing equations can be written: 


dJj _ Ome tO 1 VV, o GE, “WE ) Loa | 
az es mV, | Jy “| ee OG V2 "x | Jz 
Z 



































V V 
Pez c 2 Z 
, 2 Me | LeO Me ov, Ato X me , 2a (vx Vix) 32 
1¢vz “zo S Vz 876 2 Y 23) iat “26 
dJ 0 V, -g. ltov os k V | 
on = |e Z ren «KI J oe a | 
Z | 1°v,2 ov,¢ Oo viz «| 3 O Vz 
, IVe ‘ l+o° Vy " 0Vx arte) 7 dVz , 20 (vz-V,)91 
ave V x6 Oo Vv, IVx6 0 IVx6o eve IVxO 











Q, OV iar) 5 OV. l+o OV 20 (v_-V,)91 
=f: = eG ees Keo aoe, Z, Yd faces 
ips 5 Ge k oie O Viz " oe Oi. 














Az- —& _ ito KV, - Bez — Ginga) ee <x, 
“Gain O Vv, 1 Vz2 fo) Vz 
E. g Iason ie \ Caler So Seer 
Ce -_O& gee EIS =, aan 5 D = 
12y 2 ie Ov, 2 0 Moe 0 Vx 
: ZV 1+ 
Pee elo OV kK AV _ lite 2y7 
1y, yo Ve cee 
p2- KV, - Ito vy kV, | lto k 4V,, 
iV, O , O 


oO” 





ny Ais) 


eZ 


_— 


The six equations can 


Qu 
N} 


av 
N | 
| 


Cu 
N | 


fé 
n| 
NH 


Taking the moments of 


qy 
> 
sf& 
ia 
tt 


_ Ad2J, 


AJ, Jq + BJ35q + EJ,J5 * q, a1 Ja 


sta eal) 
5 


| 2a 


q? - (vz -Vy 


l° vz 
now be written as: 


Ain Baliae 2 ie Jig qj 


AJ, + BIg 4 EJg y 4, al 


CJz + 


CJq + 


xX 
Jy = eee 


eels 


tne whirS ter Our sedtua tiles 


d ade 
a, (3144) “es 


Cee 


eee Cyd. = 


dae) as 
3 V2.0 A yan ta5 


me eee ot Pode + qeel J, = d (J. J zi -J dJ 
: “Wyo 7 de 3a 


XO 
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MpeeecI,3,  Dd,J, . BI 
eee oe 2 PIG, ee Se tg) Fg sta 


Now, combining the above equations give the result: 


dJ | | a pale ape | oes Si imeaie® cane 
— 3. 77 8V 270 : : | —. | . 
XO 


where 


e 2 ECJgJ, -JzJ¢) i F(JjJ_ - J5d0)- 


Using the equation for 1 that was developed earlier and taking 


mie partial derivatives gives: 


O 


| De 
Bip 2 eilent( 2i(ieation) ~~ J, sa (1-z/a,) "1/3 
OV a] a 
and 
=k 5 


| aL ; Pe, 2'(1-z sfay) ST a | a- z/ 0) 


If gas diffusion had been neglected in the development, the 
above two equations would have been zero and the solution to 
Milewtransport Cquation would become unphysical. 

oct tUtinemmne sbove results imeo themda.1crencenal 


equatron for J results in: 


_(A#C-E)J + are? z'(1-z io)" 2/3 (5 a 13,42" 


a ct | 


ee 
ei: c/OMemg) | gegen eca an Sus. edee 
4 a z'(1-z ey \C 194) - J253')-(51 "54 - J2'53)) dz 
cae v, 


+Q 


S2 





where now, 


Z 
E = q, SEE 2! (1-2'/aq)°*/ aa) 


Alo v2 


Foresimpiaiication and)as a firse€ approximation assume: 


ae Ge) 
and the solution to this differential equation is: 
Z 
J = EXP | Cie oa 


O 


Seance J = 1 at Z=o. 
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